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Abstract: In February 2007 the MAGIC Air Cherenkov Telescope for gamma-ray astronomy was fully 
upgraded with an ultra fast 2 GSamples/s digitization system. Since the Cherenkov light flashes are 
very short, a fast readout can minimize the influence of the background from the light of the night sky. 
Also, the time structure of the event is an additional parameter to reduce the background from unwanted 
hadronic showers. An overview of the performance of the new system and its impact on the sensitivity of 
the MAGIC instrument will be presented. 
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Introduction 

MAGIC is an Atmospheric Cherenkov Telescope 
situated on the Canary island of La Palma 
(28.75°N, 17.86°W, 2225 m a.s.l.). Fundamental 
parameters of the telescope are a 17 m (0) mirror 
of parabolic shape (which preserves the time struc- 
ture of the Cherenkov light flash), and a hexago- 
nally shaped camera of 576 hemispherical photo- 
multiplier tubes (PMT). The field of view of the 
camera is w 3.5° 0. The fast PMT analog signals 
are routed via optical fibers to the DAQ-sytem elec- 
tronics where the signals are digitized and saved on 
disk. Further details can be found in [2]. 

Until the end of January 2007 the signals were dig- 
itized by a 300 MHz FADCs system. On February 
2007 the data acquisition was upgraded with ultra- 
fast FADCs capable to digitize at the ultra-high 
speed of 2 GSample/s [4]. In the old system, an 
additional stretching of the pulses was necessary 
in order to ensure a proper sampling of the signal. 
After the upgrade, the stretching is no longer nec- 
essary and the width of the pulses is now ~ 2.3 ns 
FWHM, half of the former value. The new sys- 
tem enhances the telescope performance basically 
for two reasons: a reduction in the amount of NSB 



(Night Sky Background) light integrated with the 
real signal (due to a smaller integration window), 
and the possibility to reconstruct with a good time 
resolution the timing characteristics of the show- 
ers. 



Signal extraction and Image Cleaning 

Currently a cubic spline is applied to find the max- 
imum pulse within the useful range of the FADCs 
(~30 ns), and the position of the maximum of the 
spline is taken as arrival time of the signal at that 
pixel. The intensity of the signal is then obtained 
by integrating the spline in a range of 7.5 ns. The 
digital filter method, used for the signal extrac- 
tion for the analysis of MAGIC data in the past, 
has not yet been tuned to the characteristics of the 
new hardware. The time resolution of each pixel 
has been estimated to be around 0.4 ns RMS for a 
40 phe signal through the study of the pixel time 
spread in calibration events, although this value 
may improve with the use of a more sophisticated 
pulse reconstruction. 

The pixels which belong to the shower image are 
selected from the whole camera picture by the so- 
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Figure 1: Cleaning features on low energy 7-MC (with a power law spectrum of index -2.6) depending on 
charge threshold (x in phe) and time constrain (y in time slices: 1 t.s. = 0.5 ns) levels set. 



called Image Cleaning (IC) algorithm. The proce- 
dure consists in setting a threshold signal value to 
select the so called "core pixels" (default 10 phe) 
and a second threshold to select or reject the neigh- 
bors as boundary pixels (default 5 phe). The idea 
behind the use of a further time constrain (impos- 
ing a limit on the arrival time difference between 
pixels) is to lower the pixel threshold, in order to 
have images with higher pixel number. In other 
words, a pixel is considered a boundary pixel if its 
charge is above the boundary threshold, at least one 
of its neighbors is a core pixel and its arrival time 
is within ±At of that of the core pixel. Requiring 
this time coincidence allows keeping more infor- 
mation from boundary pixels, avoiding to confuse 
NSB signal with real image tails (since Cherenkov 
pulses are very short in time, the probability to in- 
clude a fluctuation is very small). 

IC levels have been optimized in this study. Plots 
in figure 1 represent the characteristics of differ- 
ent levels of cleaning. Only low energy events 
have been selected, through an upper cut in the sig- 
nal of the two highest pixels (SIZE-2, a parameter 
which will not change by varying the tail cuts), 
since the main goal of the reduction of the tail 
cuts is to lower the energy threshold of the gamma 
event sample after image cleaning. The entry on 
the x axis (from 4 to 10) is the core pixel thresh- 
old imposed. The threshold for boundary pixels 
(not written) is always set as half of the core one. 
On the y axis (from 1 to 7) we have the second 
important setting level: the time coincidence win- 
dow imposed between core and boundary pixels 



(this is expressed in digitization time slices where 
1 t.s. = 0.5 ns). 

The first plot of figure 1 represents the mean num- 
ber of boundary pixels in the images. This num- 
ber increases fast by lowering the charge thresh- 
old level but a tighter time coincidence constrain 
keeps the number of boundary pixels under con- 
trol. The second histogram describes the fraction 
of the triggered low energy gammas surviving the 
IC (at least two neighbor core pixels and three pix- 
els in total are required). More events means essen- 
tially a lower energy threshold (more low energy 
events have been kept). The third plot is meant to 
classify the goodness of a certain cleaning config- 
uration. The plotted value is the naive quality fac- 
tor Q (enhancement of the signal to noise ratio) of 
the best cut in the ALPHA parameter. For a point- 
like 7-ray source, an upper cut in this parameter 
rejects a good fraction of the isotropic background 
while keeping most of the signal. In most of the ex- 
plored parameter range the quality of ALPHA does 
not change significantly. The sudden worsening of 
the image orientation reconstruction happening at a 
cleaning levels (4;2) phe is likely originated by the 
introduction of clusters of noise pixels (secondary 
islands). 

Time-related image parameters 

In order to exploit the better time resolution of the 
new FADC system, new time-related image param- 
eters have been introduced an characterized. The 
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Figure 2: TIME GRADIENT - DIST correlation for MC- 7-ray (left), real- 7-ray (center) and OFF (right). 



TIME RMS parameter estimate the arrival time 
spread of the Cherenkov photons in the pixels be- 
longing the cleaned image (//, 7-ray and hadron 
induced showers may have different characteris- 
tic time spread [10]). Previous studies [5] deter- 
mined that along the major axis of gamma induced 
Cherenkov image is present a time structure. This 
can be well approximated as a linear TIME GRA- 
DIENT. The slope of this gradient is related to 
the angle between telescope and shower axes and 
moreover, for showers of a given direction, is well 
correlated with the Impact Parameter (IP) of the 
shower. 

The possibility of using timing (multiplied by c) as 
third space coordinate of the charge distribution of 
the shower is at the moment still under investiga- 
tion. 

Monte Carlo (MC) reliability is of extreme im- 
portance to perform a good gamma/hadron sepa- 
ration and energy reconstruction. Both tasks are in 
fact done comparing simulated 7-ray to real data 
events. The 7-ray-MC database has been updated 
accordingly to the new digitization features (for ex- 
ample the value of the time jitter of the PMTs used 
in the simulation has been tuned using the cali- 
bration events). Concerning the standard Hillas[6] 
parameters the agreement between simulation and 
real events was known to be good[7], but also a 
comparison regarding the timing features is needed 
to ensure the goodness of the MC used. 

A data sample of 5.6 h of Crab Nebula observation 
(8th 10th 16th and 18th of February 2007; zenith 
angle between 5° and 30°) performed in wobble 
mode 1 have been used for this purpose. The distri- 
bution for the real data gammas is obtained as sub- 
traction between ON and OFF distribution of the 
parameter considered. A very loose cut in the 7/h 



separator parameter 2 HADRONNESS (<0.8) was 
applied to preliminarly clean the distributions from 
clear hadron candidates (while removing nearly no 
gammas) events. The agreement turns out to be 
good (see figure 3). 

A proof that the TIME GRADIENT - DIST corre- 
lation is present also in the real data 7-rays and not 
only in the MC is given in figure 2. The left plot of 
the figure is done with pure 7-ray-MC events while 
the plot on the center is the difference between 
ON and OFF distribution in the two-dimensional 
space considered. A correlation TIME GRADI- 
ENT - DIST is present in both cases whereas the 
OFF data shows no correlation (plot on the right). 
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Figure 3: Time-related image parameters, real 
data-MC comparison. A cut greater than 150 phe 
in the SIZE parameter has been applied. 



Results and Conclusions 

Three different analyses of the same previously de- 
scribed Crab Nebula data have been performed: 



1. Source is not pointed directly but 0.4° away, in or- 
der to have simultaneous ON and OFF data [3]. 

2. The Random Forest (RF) [1] 7/h separation method 
provide a single parameter to estimate the probability of 
a shower to be 7-ray or hadron initiated. 
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Crab Nebula, Size2>100, Epeak -280 GeV 



standard cleaning 10-5 
sensitivity —2.5% of Crab in 50 h 



ALPHA [deg; 



Crab Nebula, Size2>100, Epeak -280 GeV 



^ time cleaning 6-3 + time parameters 
sensitivity of Crab in 50 h 




Figure 4: Comparison of two a-plots between analysis 2 (left) and 3 (right). Keeping the same number of 
excess events background is reduced by about 50%. Energy threshold: 280 GeV. 



1. Using 10-5 IC and standard image parame- 
ters for 7/h separation (reference analysis). 

2. Using an IC of 6-3 with also the time con- 
strain (level 3 t.s.). The same standard pa- 
rameters of analysis 1 for 7/h separation. 

3. Using the same time IC of analysis 2 and 
in addition to the standard parameters the 
TIME RMS and TIME GRADIENT were 
used for 7/h separation. 

In figure 4 we compare two a-plots obtained from 
analysis 2 and 3. Time parameters allows ^50% 
better background suppression keeping the same 
amount of excess events. Similar improvements is 
seen at all energies. 

In figure 5 is shown a more general comparison: 
The telescope sensitivity 3 as function of the num- 
ber of excess events in the signal region for differ- 
ent HADRONNESS cuts. The ALPHA was opti- 
mized for each point independently. 
Improvements (of 15%) have been found also in 
the event energy reconstruction. In fact the TIME 
GRADIENT gives information about the real IP of 

I Sensitivity vs Number of excess | 



the shower and therefore it helps to distinguish dis- 
tant high energy showers from closer, low energy 
ones. 

We concentrated here on the a-analysis approach 
but further studies using the 9 2 method are already 
planned. How to apply this method in a source 
position independent framework and how timing 
could be use for the optimization of DISP are the 
first issues scheduled. 
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Figure 5: Sensitivity vs number of excess events 
(£ ; peafe=280 GeV). Cross circle and star points are 
refereed to analysis one two and three respectively. 



3. Defined as the flux that gives an excess equal to five 
time the square root of the background. 



